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Abstract 


This paper provides a detailed review of renewably driven hydrogen systems and modeling 
approaches applicable to these systems that have been reported over the last two decades. Several 
renewable energy technologies, including solar photovoltaic, wind, and hydro, have been considered 
as the power source in these reports. Storage is an important aspect of hydrogen systems, and options 
for this are summarized here. Utilization systems include fuel cells as well as a variety of thermal uses. 
Some of the reported studies have addressed residential applications whereas others were related to 
commercial scale systems. This paper particularly emphasizes aspects of modeling of the various 
components for the renewable hydrogen system. Based on the literature on this area, conclusions are 
provided on the current understanding as well as future work related to this topic. 
© 2007 Elsevier Ltd. All rights reserved. 
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1. Introduction 


As concerns grow about degradation of the environment, availability of conventional 
fuels, and the political viability of nuclear power, alternatives for energy supplies 
need to be examined. One of the hopes for the future is that renewably based approaches 
for energy can be used as replacements for our currently more-widely used sources. 
However, many of the renewable energy sources have some degree of temporal 
output, causing problems for power-on-demand requirements. For this reason, 
some form of storage will undoubtedly be of value. Also, there will be some future 
requirements for a transportable fuel that may or may not be able to be satisfied 
by batteries. 

One approach to be considered is use of hydrogen as a significant element in our energy 
future. Major advantages of using hydrogen include: firstly, it is a clean source of energy; 
secondly, it can be used for seasonal storage of energy; thirdly, it can be transported (say as 
an aircraft fuel), and finally it can be used to produce electrical as well as thermal forms of 
energy. 

The concept of integrating renewable energy with hydrogen systems was given some 
serious consideration in the 1970s [1,2]. Fischer [3] provided an interesting review on 
experimental and analytical studies of early photovoltaic/electrolyzer systems used for 
hydrogen production. Most of the experimental studies reported were on a laboratory scale. 

In last three decades, intense research has been conducted around the world on this 
topic. Emphases are not only to improve the performance of existing hydrogen production, 
storage and utilization technologies, but also to integrate them effectively with renewable 
energy sources. As noted above, one of the key elements of this is that it provides a good 
buffer to time-varying energy inputs. 

Fig. 1 shows a generic renewable hydrogen system. The renewable energy source can be 
of any type, including time-varying forms such as solar, wind, and hydroelectric. Hydrogen 
can be produced by electricity from these kinds of systems as well as by chemical and 
biological means [4]. 

Storage is an important aspect of this type of system. This is particularly true if this is 
driven by a time-varying energy source but it is also needed if there is a widely varying time 


S.S. Deshmukh, R.F. Boehm / Renewable and Sustainable Energy Reviews 12 (2008) 2301-2330 2303 


Nomenclature 

Apv area of PV module (m°) 

a modified ideality factor (a = npkTpNs/q) 

Acs be empirical constants for cell temperature calculation 

Afe active area of cell (cm?) 

age operation time for fuel cell (h) 

ref reference modified ideality factor 

Cit proton concentration at the cathode membrane/gas interface (mol/cm*) 

Ch, liquid-phase concentration of hydrogen at the anode/gas interface 
(mol/cm?) 

Cho water concentration at the cathode membrane/gas interface (mol/cm*) 

CO, concentration of oxygen at cathode membrane/gas interface (mol/cm?) 

Cow thermal capacity of cooling water (J/K) 

Cinh effective heat capacity (J/°C) 

Cy specific heat at constant pressure (J/mol K) 

Ciel overall heat capacity of electrolyzer (J/K) 

F Faraday’s number (96485) 

Gnocr instantaneous solar irradiation at SRC (W/m?) 

Gr instantaneous total solar irradiation on PV (W/m°) 

i current density draw from fuel cell (A/cm?) 

La electrolyzer current (A) 

Ir ref reference light current (A) 

I, light current (A) 

ip maximum current density (A/cm?) 

Im,ref reference maximum power point current of PV cell (A) 

Irn maximum power point current of PV cell (A) 

Io,ref reference diode reverse saturation current (A) 

Io diode reverse saturation current (A) 

Ip current of PV cells (A) 

Liere reference short circuit current (A) 

Tye short circuit current (A) 

k Boltzmann’s constant (1.381 x 107” J/K) 

ee k? reaction rate constants for anode and cathode, respectively (cm/s) 

keen empirical term to represent the apparent rate constant for anode and 
cathode 

kpr empirical term to represent the change in catalytic activity with age (V/h K) 

key effective heat loss coefficient (W/°C) 

L membrane thickness (cm) 

LMTD log mean temperature difference for heat exchanger (K) 

n number of electrons transferred per reaction (2) 

ne number of electrons involved in the reaction 

Neel number of cells in electrolyzer stack 

Ne fe number of cells in fuel cell stack 

Ngas hydrogen flow rate (mol/s) 


nH, molar flow rate of hydrogen (mol/s) 
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NO 
no, 
Npoly 
Npv 
Ns 
Mst 
P 
P a ht 
a wt 
Pe,in 
Pc,out 
P comp 
Peq 
P H2 
P 02 
P ref 
Pst 
P tho wt 
q 
Qeool 
Qgen 
Qross 
Omn 
R 


Fm 
peroton 
Ry ref 
Rs 
Rsn, ref 
Rsn 

R, el 

S pv,ref 


molar flow rate of water (mol/s) 

molar flow rate of oxygen (mol/s) 

polytropic index 

parameter 

number of cells in series 

number of moles of hydrogen inside the storage tank (mol) 
operating pressure of electrolyzer cell (N/m?) 

actual available power of hydro turbine (W) 

actual available power of wind turbine (W) 

pressure at the outlet of compressor (N/m*) 

pressure at the outlet of compressor (N/m*) 

power required by compressor (W) 
equilibrium pressure of hydrogen gas (N/m’) 

partial pressure of hydrogen (N/m?) 

partial pressure of oxygen (N/m”) 

reference pressure (N/m7) 

pressure inside the storage tank (N/m?) 

theoretical available power of wind turbine (W) 

electronic charge (1.609 x 107° As) 

heat lost to cooling media (W) 

heat generated by elecrolyzer (W) 

heat loss (convection and radiation to ambient) by electrolyzer (W) 
total heat production/consumption (W) 

gas constant 


ri, F2, S, ty, to, t3 empirical constants for electrolyzer 
Releetronic 


resistance to flow of electrons at the electrodes (Q) 
membrane specific resistivity for the flow of hydrated protons (Q cm) 
resistance to the flow of protons in the cell (Q) 
reference series resistance (Q) 

series resistance (Q) 

reference shunt resistance (Q) 

shunt resistance (Q) 

overall thermal resistance of electrolyzer (W/K) 
reference effective absorbed solar irradiation (W/m?) 
effective absorbed solar irradiation (W/m?) 
operating cell temperature (K) 

ambient temperature at SRC (K) 

ambient temperature (K) 

cooling water temperature at inlet (K) 

cooling water temperature at outlet (K) 

back surface temperature of PV module (K) 

cell temperature at standard rated condition (SRC) (K) 
cell temperature of PV (K) 

reference temperature (K) 

standard cell potential (V) 

standard cell potential for fuel cell (V) 
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Unx 


Uai 
Ure 
UŁ nocr 


Nespv 
"comp 
"conc 


Nohm 


a, B, y, 6, e empirical constants for calculation of specific heat at constant pressure 


Či, Č2, 3, €4 empirical parameters for fuel cell 


product of overall heat transfer coefficient and area of for heat exchanger 
(W/K) 

cell potential for electrolyzer (V) 

overall cell potential for fuel cell (V) 

overall heat transfer (convection and radiation) coefficient of PV at SRC 
(W/m? K) 

overall heat transfer (convection and radiation) coefficient of PV (W/m? K) 
reference maximum power point voltage of PV cell (V) 

maximum power point voltage of PV cell (V) 

reference open-circuit voltage (V) 

open-circuit voltage (V) 

voltage of PV cells (V) 

reversible cell potential (V) 

thermo-neutral potential (V) 

volumetric capacity of storage tank (m°) 

wind speed (m/s) 


chemical activity parameter for the cathode 

change in Gibbs free energy (J/mol) 

standard state free energy of the cathode reaction (J/mol) 

standard state free energy of chemisorption from the gas state (J/mol) 
change in enthalpy of formation (J/mol) 

standard (10° N/m? and 298 K) enthalpy of formation (J/mol) 

change in entropy of formation (J/mol K) 

standard (10° N/m? and 298 K) entropy of formation (J/mol K) 
temperature difference between cell and module back surface at a 1000 
W/m? irradiance (°C) 

activation losses (V) 

module efficiency to convert incident solar irradiation to electrical energy 
compressor efficiency 

concentration losses (V) 

ohmic losses (V) 


variation of use. This can be accomplished using compressed gas storage tanks, liquefied 
hydrogen storage, metal hydrides or carbon-based materials. We also consider short-term 
electrical storage as an option. This could be batteries or grid connection, if the latter is available. 


Renewable 
Electrical Source 
Hydrogen p ER Short Term 
Generation,Storage 4 Storage 


Load (Thermal l4 


and Electrical) 


Fig. 1. Generic renewable hydrogen system. 
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Loads for houses can be classified on an electrical or thermal basis. Electrical loads can 
be satisfied by fuel cells (stationary and automobile application) whereas hydrogen can be 
combusted for satisfying thermal loads, e.g., IC engines (cars), cooking, water heating, 
HVAC systems, etc. 

A variety of studies that have been performed on combined renewable energy and 
hydrogen systems can broadly be classified as follows. 


1.1. Renewable hydrogen generation system (RHGS) 


The system generates hydrogen from renewable energy sources and can be either 
stand alone or grid connected. This system is a typical hydrogen refueling 
station. Note that for residential applications such a system can fulfill the hydrogen 
requirements for fuel purposes only. A description has been given for the initial 
design, operation and cost optimization of the system with photovoltaic cells connected 
to an electrolyzer [5-9]. Experimental data on small-scale systems with proton 
exchange membrane (PEM) electrolyzers have been summarized [10]. Experimental 
data for RHGS with PV that was installed on a house in Switzerland have been 
reported [11]. 

One of the first studies including description of component models to simulate the 
performance of the RHGS was given by Carpetis [12]. Others have described simulations 
performed on the TRNSYS software for similar systems [13]. Hydrogen production as a 
function of maximum power point and tilt angle has been given. The performance 
modeling of the RHGS (with PV) was reported [14]. Further, cost optimization studies of a 
photovoltaic electrolyzer system using HYSOLAR along with performance simulations 
were reported by Seigel and Schott [15]. Other cost optimization studies have also 
appeared in the literature [16,17]. 

Some of the fairly recent studies discussed the preliminary system design and economic 
feasibility of a combined PV and PEM electrolyzer [18,19]. Bilgen developed models to 
study performance and cost optimization of PV electrolyzer systems [20]. The paper 
reported comparisons for fixed and tracking panels. Others have reported the optimum 
hydrogen cost for a small-scale (< 100 W) directly connected RHGS with PV [21]. Another 
study has focused on the use of PSCAD/EMTDC software for the design and simulation 
of transient performance of RHGS with PV [22]. 


1.2. Renewable hydrogen generation and utilization system (RHGUS) 


This system not only generates hydrogen from renewable energy sources, but also 
includes a load that utilizes hydrogen. The load this system serves could be electrical and/ 
or thermal in nature. The system can be either stand alone or grid connected. Such a 
system includes on-site/off-site hydrogen generation also. The latter refers to a system in 
which hydrogen is generated at a renewable source facility and transported by some means 
to the point of use. 

Design considerations based on efficiencies and energy calculations for various 
components of RHGUS with PV as an energy source were reported [23,24]. The long- 
term operating experience for large-scale demonstration plants of RHGUS systems 
installed in Germany has been described [25,26]. Griebhaber summarized the energy 
consumed and generated for RHGUS installed on self-sufficient Solar House [27]. The 
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efficiencies for various components of RHGUS that supplied continuous power to an air 
compressor have been noted [28]. Experimental data for the performance of a combined 
PV and wind turbine-operated RHGUS for remote applications were summarized. A 
programmable load was used to provide the energy demand for the house [29]. Others 
reported the performance of RHGUS with PV installed to supply continuous power 
(200 W) for a remote application. A metal hydride bed was used to store the hydrogen 
generated by an electrolyzer [30]. 

Models were developed and overall system simulations were performed to optimize the 
size of PV array area for the RHGUS system [31]. The program H2PHOTO was 
developed and simulations were performed based on simple models (electrical and 
thermal) to design and predict the performance of 1.3kW RHGUS system [32]. TRNSYS 
was used to simulate the annual performance of RHGUS system for residential 
purposes. A catalytic burner was used for providing heat to the house but the sizing of 
the components depends on both the electrical and thermal loads of the house [33]. 
Vosen and Keller [34] described time dependent annual simulations for the RHGUS with 
PV for stand-alone applications using state of charge and neural network control 
systems. The cost comparison between hydrogen only, battery only and combined 
system was also provided. Others developed a model to estimate hourly solar radiation on 
a tilted surface from monthly average daily insolation. This was further used to 
estimate the energy available from a PV (model based on efficiencies) and wind turbine 
combined power source for a RHGUS [35]. The RHGUS system for stand-alone 
residential application was simulated for PV only, PV-—micro-hydro turbine and 
combined PV, wind and micro-hydro turbine as energy sources [36,37]. A thermo- 
economic analysis was performed for RHGUS with PV [38]. The optimization to 
minimize the investment cost of RHGUS with PV and micro-hydro turbine was carried 
out [39]. The results gave local minima but it was not sure that the model gave global 
minima. Short-term performance of the RHGUS with PV and wind turbine for the 
control system and power conditioning devices developed by authors was also reported 
[40]. For the similar system, the effect of operation of electrolyzer with oxygen recovery 
on the performance of overall system has been described [41]. The dynamic models of 
various components of RHGUS with PV and wind turbine were developed and 
compared with experimental results. Further, these models were used to simulate 
the system performance for residential purposes [42]. The importance of the control 
strategies, based on simulations performed on TRNSYS, to operate the RHGUS 
with PV has been noted [43]. Annual simulation for RHGUS with PV and wind 
turbine was performed for an average load of 1kW using Hybrid2 software developed 
by National Renewable Energy Laboratory (NREL) [44]. Performance simulation 
of the RHGUS system with PV and reformer, installed in Stockholm, Sweden, to 
obtain hydrogen from biogas has been outlined in the literature [45]. Shapiro et al. 
reported the performance of an electrolyzer (up to 200 psig) placed inside the storage tank 
that can be combined with RHGUS [46]. Dynamic modeling of hydrogen side 
(electrolyzer, storage, and fuel cell) of RHGUS has been performed for application for 
residential application when combined with PV array [47]. The annual performance 
modeling of RHGUS with PV for a standard residential house in Las Vegas, NV has been 
performed by Deshmukh and Boehm [48]. The hydrogen needed for fuelling the 
automobile and other residential usage (water heating, clothes drying, cooking, etc.) has 
been considered. 
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2. Modeling by component 


This section reviews various approaches to modeling of various aspects of renewable 
hydrogen systems. Models considered are based upon their applicability for system 
simulation purposes. Both components and input phenomena are included here. 


2.1. Input resources and hourly data 


Modeling of various components requires a variety of input data. Hourly averaged 
annual data for solar irradiation, temperature, wind speed, wind direction, and water flow 
are generally needed for modeling. 

The hourly averaged solar irradiation data for US can be found on Renewable Resource 
Data Center developed by NREL [49]. Wind data can be downloaded from the NREL 
website [50]. Monthly wind rose chart information averaged over thirty years can be 
obtained from United States Department of Agriculture website [51]. This can be helpful 
for the design of wind turbine applications for a particular location. For annual 
performance, hourly wind speed and direction data are needed for a particular location. 
Modeling of hydro turbines needs hydrograph data which are listed on the United States 
Geological Survey website [52]. This data can help to design the micro-hydro turbines and 
predict their annual performance. 

Hourly data obtained can be further used to calculate the amount of solar irradiation 
incident on the PV panels. The amount of incident irradiation also depends on the type of 
mounting system used for PV panels. Sun-tracking systems generally have higher incident 
energy compared with that on a fixed system planar system. Tracking systems can be 
classified as follows: 


è Single axis tracking: 
O horizontal east-west, 
O horizontal north-south, 
O vertical. 

è Two-axis tracking. 


The amount of energy available to a variety of collector types and orientation has been 
developed for 239 physical locations in the US by the NREL [53]. In doing this, hourly 
data averaged over 30 years were used. Hourly data are available for simulation purposes. 
One comprehensive source is the TMY2 data [49]. 

Sun-tracking systems can enhance the collection yield significantly. When a large 
portion of the radiation is direct (compared with diffuse), concentrating systems can be 
used. The relation to calculate the amount on incident solar irradiation based on type of 
collector, whether tracking or not, is explained in a number of sources including the 
excellent text by Duffie and Beckman [54]. 


2.2. PV modeling 
Photovoltaic cells are currently available in single and multi-junction forms. A problem 


inherent in photovoltaic devices is that they are able to interact only with a portion of the 
solar radiation distribution. Multi-junction cells are able to capture a wider portion of the 
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spectrum, and they show higher conversion efficiencies from sunlight-to-electricity than do 
single junction cells. Commonly available photovoltaic panels use single crystalline, 
polycrystalline, thin film or multi-junction cells. Many researchers have studied renewable 
hydrogen systems with photovoltaic cells (as mentioned in previous sections) as a source of 
energy for generating hydrogen. 

Modeling PV arrays for renewable hydrogen systems can be accomplished by using the 
single diode model with equivalent circuit given in Fig. 2 [54]. This model uses the data 
provided by manufacturers (Ise Uso Imps Ump Hisee and Huoc) to model the PV 
performance. 

The IV curve is calculated as a function of incident solar irradiation and cell temperature 
using Eq. (1) as follows: 


Ip = 11 ~ Lo exo | 
a 


Upo + TpoRs 
Rsh 


Up T rE] oo | 


(1) 


where light current (Iz), diode reverse saturation current (Io), series resistance (Rs), 
shunt resistant (R,,), and the parameter a are the five characteristics that are needed to 
solve Eq. (1). Based on the data from manufacturers, Eqs. (2)-(7) are solved 
simultaneously to obtain the reference values: 

At the short circuit condition, J,, = Ise and Up, = 0, 


Tserep Rsre 
Lseref = Th yref = LO,ref G Pa a | 


Aref 
z I seref Rs ref (2) 
Rshyef 
At the open-circuit condition, Zp, = 0 and Up = Uoc 
U Uden 
0= Let = Iosef lexp | set | oc,ref (3) 
dref Rshref 


l l 
Oi wt oe f 


Fig. 2. The single diode equivalent circuit model. 
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At the maximum power point, 


Um + m Rs ref 1 
dref 


In = Thref = Loref lexp | 
_ Um + TnRs ref f 


z (4) 
sh,ref 
Also at maximum power point, 
dP pv dI pu 
dU = 0=Im + Uny i (5) 
P’\lm PY\lm 
where 
Ioref exp U mref tLm ref Rsref 1 
dI pv ref dref Rshref ( 6) 
dU = Lose Rsref U mref +L mrep Reset Reses ` 
po ef ss z Z 3 
K 1+ Aref exp ref + Rshref 


And, checking the condition for estimation of the open-circuit voltage predicted by 
model, it is given as 

OU ce =u wy Uo T pv) ae Ud T pe ref) 

OT al To= Tre 


(7) 


The reference values are used to calculate the five parameters at operating condition 
using Eqs. (8)-(13): 


Tye 
pea E, (8) 
Pr ref 
Spo 
Ir = bee + hr (T pv = T pog), (9) 
pv,ref 
T. i E E 
AE = a E = tell (10) 
f (z=) kT oe kT T. 
Here, 
Eg = Egl! — CT p — Tpv,ref)l (11) 
S vre, 
Ror = Rory L; (12) 
pv 
Rs = Rs, (13) 


Finally, using the 5 parameters at the operating conditions, we can either calculate the 
current at a given voltage or use (14) and (15) for calculating current and voltage at 
maximum power point: 


m ImRs 
yD pee E lexp Ante’ = | 
a 


Um + ImRs 


Rsh ( : a 
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2exp [Eatin Re] Ri 


ee Tak exp [Eatin Re ae n 


0=In+ Um (15) 


sh. 


Performance of the PV cell is very much dependent on the cell temperature which itself is 
a function of solar irradiation, ambient temperature, wind speed, wind direction type of 
mounting, and other factors. An energy balance for the PV cells can be used to find the cell 
temperature from Eq. (16) [54]: 


T po m Ta = Gr ULNocr i ten] (16) 
Tyoctr —~Tanocr Gynocr UL (tor) |” 
where 
In Um 
= . 17 
Inpe ApyGr ( 


Eq. (16) can be solved in an iterative way to obtain the cell temperature. The cell 
temperature can also be solved as an explicit function of incident solar irradiation, ambient 
temperature, and wind speed using Eqs. (18) and (19) [55]: 


Gr 


Tv = Tmp +>— 
a P? Gyocr 


AT (18) 


where T, is the back surface temperature of module and is given by 
Tipi =T, + Grethe WS). (19) 


Figs. 3—6 show the performance of PV panels of 65 W rating. The data available from 
the manufacturer are e= 4.5V, U,.= 21.4, Im = 3.95, Um16.5, Hīse = 0.00026, 
Huoc = —9.85, Ap, = 0.633, and Ns = 36, as reported in [54]. Using the 5-parameter 
model noted above, the performance of the PV panels has been calculated. Fig. 3 indicates 
the -U curves at different levels of solar irradiation at a cell temperature 40°C. The 
maximum power point curve is also shown in Fig. 3. -U curve and the power vs. U curves 


Current (Amps) 


0 5 10 15 20 
Voltage (Volts) 


Fig. 3. I-U curve and maximum power point curve for different incident solar irradiation (W/m?) at 40°C for a 
PV panel. 
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5 70 
7 60 
= 7 50 _ 
Q £ 
£ J 40 § 
z pa 
2 —I-VCurve 7 30 g 
3 a 
oO Power-VCurve + 20 
7 10 
0 1 1 1 1 0 


0 5 10 15 20 25 
Voltage (Volts) 


Fig. 4. I- U curve and power vs. U curve for a PV panel at 1000 W/m? solar irradiation and 25°C temperature. 


—1000 W/m2 50 C 
—— 1000 W/m2 25 C 


Current (Amps) 
N 
a 


0 1 1 li 
0 5 10 15 20 25 


Voltage (Volts) 


Fig. 5. The effect of temperature on the performance of a PV panel. 


50 L — 1000 wim2 
——500 Wim2 


Power (Watts) 


20 25 30 35 40 45 50 
Temperature 


Fig. 6. Effect of temperature on maximum power of a PV panel. 
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Table 1 

Values of empirical constants a,, be, and AT,, for various types of PV panels and mounting arrangements [55] 
Module parameters Module Ac be AT,, CO) 
Glass/cell/glass Open rack —3.47 —0.0594 3 
Glass/cell/glass Close roof mount —2.98 —0.0471 1 
Glass/cell/polymer sheet Open rack —3.56 —0.075 3 
Glass/cell/polymer sheet Insulated back —2.81 —0.0455 0 
Polymer/thin film/steel Open rack —3.58 —0.113 3 

22X linear concentrator Tracker —3.23 —0.13 13 


are shown in Fig. 4 at a solar irradiation of 1000 W/m? and cell temperature of 25°C. The 
effect of temperature on the /-U characteristics of the PV panel is depicted in Fig. 5. The 
effect of temperature on the maximum power point for different levels of solar irradiation 
is shown in Fig. 6. It can be observed that performance changes significantly as the cell 
temperature changes. It is crucial to develop the methods to keep the cell temperature as 
low as possible to improve the performance of the PV cells (Table 1). 

Other formulations for modeling photovoltaic devices have been reported. For example, 
King [55,56] and Luft et al. [57] have given examples of this. King et al. [55,56] developed 
an empirical model for simulating the performance of various PV technologies (the models 
take electrical, thermal, solar spectral and optical effects into consideration). 

Empirical parameters found were specific for each type of module and a list of these is 
maintained at Sandia National Laboratories website [58]. Luft et al. [57] developed a 
model for predicting the current produced based on the Ise, Im, Um, Uoc- Hadj et al. [59] 
developed different relations to evaluate the five parameters used in Eq. (1) for modeling 
PV performance. Comparisons of these models for representing PV performance are given 
in [60]. 


2.3. Wind machine modeling 


Studies of combined PV and wind power have been reported [26,29,35,36,40—42,44]. 
Modeling was generally based on the machine performance data available from the 
manufacturer. Fig. 7 shows the experimental and calculated (using (21)) wind speed vs. 
power generated for a 2.5kW rated wind turbine [61]. With the wind data for a location 
known, the amount of power generated can be estimated over a year’s period. 

Theoretical available wind power for a particular location can be calculated using 

Pango = PairAin WS? (20) 


wt 


The power in the wind must be modified by the “fractional efficiency” to determine the 
power out of the wind machine. Typical values of this parameter will be in the range 
0.15-0.30. 

A polynomial fit can be obtained to obtain power as a function of wind speed for any 
location and a particular wind turbine. The result for the data shown in Fig. 7 is as given in 
Eq. (21): 


Pawi = —0.177 WS4 + 3.164 WS? + 6.757 WS? — 25.042 WS + 2.119. (21) 
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Fig. 7. An example trace of power output for a range of wind speeds is shown (after [61]). 
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Fig. 8. Maximum power generated for a particular hydro turbine is shown as a function of water flow rate and net 
head [62]. 


2.4. Micro-hydro turbine modeling 


Studies of micro-hydro turbine powering hydrogen generation have been 
reported [36,37,39]. The modeling was done based on the data available from 
the manufacturers of the equipment used. Generally for residential applications, 
low head hydro turbines (e.g., reaction turbines) can be used. An example of data 
available from a manufacturer for a reaction turbine with rated power 2.2 kW [62] is shown 
in Fig. 8. 

Actual power generated for the hydro-turbine as a function of head and flow rate can be 
calculated using Eq. (22). Note that power generated is a linear function of flow rate and 
head available. The efficiency for the turbine given is approximately 


Pam = npgQh. (22) 
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2.5. Electrolyzer modeling 


Electrolyzers are devices that use electricity to split water into hydrogen and oxygen. The 
semi-empirical model of an alkaline electrolyzer was reported in [63,64]. The reactions at 
the anode and cathode of the electrolyzer cell are given as follows: 


Anode : H20 = 2H* + 50) + 2e”. (23) 


Cathode: 2Ht+ + 2e7 = Hb. (24) 


The semi-empirical model described in [63] could be easily applicable to PEM 
electrolyzers with a change in the empirical constants. The voltage required for the 
electrolyzer is obtained from the summation of the reversible voltage and summation of 
losses due to various resistances in the cell. The losses can be classified as activation losses, 
which quantify for slow electro-kinetics, ohmic losses in the cell and the concentration 
losses, which occur at higher current densities due to improper mass transfer of the species: 


Uel = Urev T Nact T Nohm Eg Neone: (25) 


Reversible voltage as a function of temperature and pressure can be given from the 
Nernst equation [64] as follows: 


U, = = [Uo + RT In [onpa] À (26) 


where Up is the standard cell potential and is calculated from the standard Gibbs free 
energy using Eq. (27). For the case of an electrolyzer reaction, the Gibbs free energy 
change is a positive quantity, and at standard conditions (25°C and 1 bar) it is equal to 
237.13 kJ/mol: 


AG 
~ nF’ 

The Gibbs free energy is a function of temperature and pressure. It can be obtained from 
a change in the enthalpy of formation and a change in entropy using Eq. (28). The change 
in enthalpy of formation and entropy as functions of temperature and pressure can be first 
obtained for individual species from Eqs. (31) and (32) [65]. Then the other parameters can 
be evaluated for the electrolyzer cell using Eqs. (29) and (30): 


Uo (27) 


AG = AH — TAS, (28) 
AH = AHy, + AH, — AH p;o, (29) 
AS = AS, + SAS 0, — ASn0, (30) 
T 
AH = AH? + f Cp dT, (31) 
Tref 
T 
as=as?+ | p= Ria | 2). (32) 
Tref T Pref 


The specific heat at constant pressure can be expressed as a function of temperature 
using Eq. (33) [66]. Constants «, P, y, 6, and ¢ can be obtained from Table 2. This function 
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Table 2 

Constants for calculating specific heat at constant pressure, Eq. (33) [66] 

Gas a Bx 10° yx 10° ôx 10° ex 10!" 
H20 4.070 —1.108 4.152 —2.964 0.807 
O2 3.626 — 1.878 7.055 —6.764 2.156 
Hp 3.057 2.677 —5.810 5.521 —1.812 


can be integrated and further used to evaluate changes in enthalpy and entropy of 
individual species using Eqs. (34) and (35): 


C, = Ria + BT + yT? + 6T? + eT], (33) 


T? yT? 6T* T’ 
ma +T), (34) 


= o 
AH = AH +R(ar4 7 + 3 re 4 5 


yT? ôT? eT’ 
AS=As?+ R(ain T+ pT +" +5 +7) 


-ra (2), (35) 
P ref 


The reversible potential U,., can be easily determined from Eq. (26) as a function of 
temperature and pressure. Following that, the electrolyzer cell potential can be evaluated 
from Eq. (36). Note that the empirical parameters r, s, and ¢ can be obtained by curve- 
fitting manufacturers’ data. The second term on the right-hand side quantifies the ohmic 
resistance: 


r t 
Uel = U rev F q. +s log (ste a ij: Gq) 


where the empirical parameters r and ¢ can be defined as function of operating temperature 
using Eqs. (37) and (38). Note that the temperature is in Celsius units for calculating r 
and t: 


ease oe (37) 
Tyre j T, j 
t=ħh+h= +% a (38) 
The amount of hydrogen produced is given by 
r Ncellel 
NH, = NF ar (39) 


where np is the Faraday efficiency which is the ratio of actual hydrogen produced to 
theoretical maximum possible production. Similarly, the amount of water required for a 
specific amount of hydrogen production and the corresponding amount of oxygen 


S.S. Deshmukh, R.F. Boehm / Renewable and Sustainable Energy Reviews 12 (2008) 2301-2330 2317 
produced are given by 


H0 = NH, = 2No,. (40) 


Performance of the electrolyzer is dependent on the operating temperature. The 
temperature can be either assumed constant or can be calculated by using a thermal model. 
A lumped thermal capacitance model [63] can be used to predict the temperature of the 
electrolyzer: 


dT 
OT = Oen _ Qinss E Qooor (41) 


where Qgen is the internal heat generation of the electrolyzer, Qjoss the heat lost to the 
ambient, and Qoo; the amount of cooling required to keep the electrolyzer at a rated 
condition. These are given as follows: 


Cie = Neel Ue — U mMer, (42) 
where 
AH 
Um = nF? (43) 
a Ce (44) 
loss — Ria a). 


A common method of cooling an electrolyzer is to flow excess water through the unit as 
compared with the water required for the actual hydrogen consumption. In this case, it can 
be assumed that the water flow related to the cooling requirement is 


Osa) = CewlT ewi E Tewo) = UA hx LMTD. (45) 


It is worth mentioning that there is a paucity of publications related to studies of 
electrolyzers. More data are required to model the electrolyzer properly based on the 
physical behavior. It is imperative to provide a good basis for the semi-empirical 
parameters used in the model for calculating I- V curve, thermal model, etc. 

Using the above approach, the effect of temperature and pressure on the /-V curves of 
the PEM electrolyzer is plotted. In the report by Cropley [67], experimental data were 
available for only one temperature and one pressure. Using this data, the empirical 
constants were obtained for the best fit and the effect for different temperature and 
pressures were calculated. See Figs. 9 and 10. 

Some recent studies have focused on the search for improved and cheaper catalyst 
alternatives for hydrogen and oxygen evolution. The performance of cobalt and nickel 
glyomines is found to be inferior as compared with platinum for hydrogen evolution at 
cathode [68]. Authors reported the performance of Ir/Sn, Ir/Ru-based oxides for oxygen 
evolution at the anode. An efficiency of 94.4% for a single cell PEM electrolyzer was 
achieved with Ir/Ru at anode and Pt/C catalyst at cathode. Ir/Sn was not found to be 
effective catalyst for oxygen evolution [69,70]. 
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Fig. 9. Effect of operating pressure on the performance of an electrolyzer cell. 
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Fig. 10. Effect of temperature on the performance of an electrolyzer cell. 


2.6. Fuel cell modeling 


A fuel cell is an electrochemical device that converts hydrogen and oxygen gas into 
electrical current (DC). A review of fuel cell modeling based on analytical, semi-empirical, 
and mechanistic approaches has been given by Cheddie and Munroe [71]. Modeling of fuel 
cells for system analysis purposes was initially proposed by Vanhanen et al. [32], 
Chamberlin [72], Lee [73], and Amphett [74]. Further improvement in the fuel cell model of 
Amphett et al. was reported in [75-79]. 

The reactions at the anode and cathode are given as follows: 


Anode: Hy = 2H* + 2e7. (46) 


Cathode : 2H* + 402 + 2e7 = H20. (47) 
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Electrical energy, water, and heat are the three products of this reaction. The total cell 
potential can be obtained from 


Ufe = U rev F Nact a "conc T Nohm: (48) 


where U,,, is the reversible voltage as a function of temperature and pressure can be 
obtained from the Nernst equation [65]. The first term gives the variation of standard 
potential of the fuel cell with respect to temperature and pressure: 


a n [Uog +RT In (papi); (49) 


Uo fe is the standard cell potential of the cell. It can be calculated as a function of 
temperature and pressure from the relations developed for electrolyzers. The only 
difference is that the products and reactants are opposite from that of an electrolyzer and 
are given by Eqs. (50) and (51): 


AH = AH n,o — AH p, — HAH o,, (50) 


AS = ASy,0 — AS, —1ASo,. (51) 


Activation polarization 77,¢; is the extra voltage applied to counteract the slow reaction 
rate at the electrodes of fuel cells. The reaction rate can be increased either by increasing 
the process gas flow rate and thereby increasing the process gas kinetic energy; or by 
adding extra voltage and thereby reducing the energy output of the fuel cell. According to 
Faraday’s law, the reaction rate is proportional to the current density; and a lower reaction 
rate implies lower current densities. Other means of reducing the activation potential 
besides increasing the process gas flow rates include: increasing active area of the electrode, 
increasing operating temperature and increasing the utilization of catalyst [73]. 

The expression for the activation polarization (anode + cathode) was derived in [78,79] 
from Butler-Volmer equation as follows: 


Nact = 1 + T+ &T In(cG,) + &4T In (i), (52) 
where či, č2, €3, €4 are given as follows: 
eo s 
b= ln [FREA Cok | 
j (Š + =) In [Aye] + (=) In [aF |; (54) 
= pla) (55) 


j R R 
&=- (5 =+ =) (56) 
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Numerical values for the €), č2, €3, €4 parameters were provided for a Ballard fuel cell 
[77,78] as follows: 


é = — 0.948(+0.004), 
E= kpk m + k + 0.000197 In (4f) + 4.3 x 10° In(c%,), 


€; = 6.8 + 0.02 x 1075, 
é, = — 1.97 +0.05 x 1074. 


These values could be changed within given limits for modeling PEM fuel cells from 
other manufacturers. 

Ohmic polarization, ohm, is the extra voltage needed to overcome the electrical 
resistance offered by the cell to the electron and proton transfer within the cell. Note that 
the resistance to the electron transfer could be due to electrodes, graphite collector plates, 
and other fixtures that connect the membrane electrode assemblies (MEA) together. This 
resistance could be decreased by using materials that are good electrical conductors. 
Proton conductivity could be improved by reducing the thickness of the membrane. The 
ohmic polarization can be estimated as follows [78]: 


Nokm = — i( Relectronic te Rees | Rinternal . (57) 


Here, R°““"°"”" includes resistance due to contact resistance. For pure graphite plates the 


actual Ronie should not be significant in comparison to R°” [79]. R’’?°" can be 
obtained as follows: 
TnL 
preron = m , (58) 
Afe 


where L and Ap are the thickness and the area of the membrane, respectively, and r,,, is the 
membrane specific resistivity for the proton transfer from the membrane. The latter is a 
function of the type and characteristics of the membrane, water content of membrane, 
temperature, and current density. For Nafion membranes the following relationship for r, 
was proposed [78]: 


181.61 + 0.03 (4) + ooe (4) | 


a [2° + Aprage — 0.634 — 3()] exp (4.18 5 323]) 


(59) 


The voltage degradation effect on the performance of fuel cells was studied in [79]. 
Aging phenomena are introduced in the model through the / terms in ohmic losses and keen 
in activation potential losses. 

Concentration polarization yon. 1s the extra voltage required to overcome the resistance 
to the mass transport of the process gases at the reaction sites on membrane electrode 
assembly. Mass transport includes the transport of gases to the reaction site (flow through 
channels and diffusion in the electrodes) and removal of water and other impurities from 
the reaction site. At the lower current densities, these losses are not significant as compared 
with higher current densities. At higher current densities, process gases should reach the 
reaction sites at higher rates and similarly the water and other impurities should be 
removed more quickly from the reaction sites. Concentration polarization can be estimated 
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RT i 
one = —a |1- >|. 60 
"leone = pF | ip ey 


Hence, the /-U curve of the PEM fuel cell can be obtained using Eq. (48) considering the 
losses due to activation, ohmic and concentration polarizations. Figs. 11 and 12 show the 
performance of the 1.4kW fuel cell stack. The experimental data used for this purpose 
were reported on the web [81]. Figs. 11 and 12 show the effects of pressure and temperature 
on the performance of the fuel cell, respectively. 

Similar to electrolyzers, a fuel cell’s performance is also temperature dependent. In a 
manner similar to what is done relative to an electrolyzer a thermal model can be used for 
predicting the temperature of the fuel cell. Internal heat generation for this case would be 
as follows: 


O rs = Ne fel Um = UfM. (61) 
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In addition, a fuel cell’s heat loss can be calculated using a heat balance across 
the inlet and outlet streams of hydrogen and oxygen flow. In the case of a hydrogen/air fuel 
cell, it must be remembered that the oxygen side balance involves 21% of the total 
atmospheric air flow. Some recent studies have focused on accurately predicting the 
activation and concentration polarization [82-84]. The evaluation of activation and 
concentration polarizations at the anode and the cathode is based on a function of the 
concentration and partial pressures of the reactant gases (H> and O2) that are related to 
each other by Henry’s law. Once the exchange current density is obtained, the 
approximations of Butler-Volmer relationship for lower and higher polarization can be 
used to obtain the activation potential for the anode and the cathode [83]. A literature 
review regarding the mechanisms of hydrogen oxidation reaction on different 
platinum electrodes and their effect on the exchange current density has been given by 
Mann et al. [84]. 


2.7. Storage modeling 


Hydrogen storage is one of the important aspects of the path towards a hydrogen 
economy. As was noted in the early part of this paper, many of the renewable sources have 
energy outputs that are quite variable in time, so a means of storing the results of capturing 
the energy is critical. The following are the commonly used hydrogen storage technologies 
[85-88] currently under research: 


1. Compressed hydrogen. 

2. Liquid hydrogen. 

3. Metal hydrides. 

4. Carbon-based materials (fullerenes, carbon nanotubes, activated carbons). 


Discussions for all these forms and their current status have been given [85,87]. In this 
section, we will focus on modeling aspects of these technologies. 


2.7.1. Compressed hydrogen storage 

Modeling of compressed hydrogen gas storage can be accomplished by assuming with 
little error that hydrogen is an ideal gas. The amount of work required to compress 
hydrogen to a required pressure can be obtained by considering a polytropic process for 
compressor. The actual work required can then be found using a compressor efficiency, 
information usually furnished by the compressor manufacturer, using Eq. (62): 


(poty—1)/Npoty 
ie Ga | (62) 


And pressure inside the storage tank can be obtained from 


Ag RT 
Pa =p (63) 


1 Npoly 


Preomp = Ngas RT 


"comp Npoly -1 


2.7.2. Liquefied hydrogen 
In general, it is not anticipated that this approach will be used for most stand-alone 
hydrogen systems. Although it is possible, it seems that the nature of a liquefaction 
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Fig. 13. Shown are the adsorption curve (upper curve) and desorption curve (lower curve) for a metal hydride 
hydrogen storage [11]. 


approach may be too complicated. Of course, the storage energy density is relatively high 
using a liquid compared with using even a very high-pressure gas. Good sources of 
information on this general topic are given in the literature [85,87,88]. 


2.7.3. Metal hydrides 

Many metals and alloys have a tendency to combine with hydrogen to form hydrides 
[89]. When hydrogen is stored inside a hydride (also called as charging or hydriding), the 
heat is released during this filling operation. On the contrary, heat input is needed while 
taking hydrogen out from the hydride, also called discharging or dehydriding. This 
reversible process for metal M can be expressed as 


M + 0.5xH <= MH, + heat. (64) 


Including a metal hydride with the renewable hydrogen system needs a fairly simple 
model that includes the rate of hydrogen storage (includes hydriding and dehydriding), the 
amount of heat generated/required, as well as the operating temperature and pressure. 
Detailed modeling of the metal hydrides is reviewed in [85]. A database of metal hydrides 
and their properties is available on the web [90]. Using the lumped capacitance method the 
temperature of the bed as a function of time can be calculated [91] as well as the heat 
generated or required for hydriding or dehydriding. Based of bed temperature, the 
operating pressure can be found using Van Hoffs equation [85]: 


Cm O- Ont) — Keg TO) — Td, (65) 
AH AS 
In eq) = pr R’ (66) 


where Q,,,, is a function of time because the heat would be released during the absorption and 
required during desorption. A basic assumption usually made in the modeling of hydride storage 
is that the cycling efficiency is 1, which means that there is no change 
in the adsorption and desorption properties of the MH. Fig. 13 shows the hydrogen adsorption 
and desorption curve at 5.3°C for a MH bed of 19 N m° rated hydrogen storage capacity [11]. 
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Metal hydride beds for hydrogen storage integrated with RHGUS have been reported in 
the literature [11,30,34]. Experimental data for the metal hydride bed integrated with an 
electrolyzer and a fuel cell were given by Vanhanen et al. [92]. 


2.7.4. Carbon-based materials for hydrogen storage 

Carbon-based materials are capable of storing hydrogen. A review of work related to 
hydrogen storage in fullerenes, carbon nanotubes and activated carbons has appeared in 
the literature [86,93]. A wide range of hydrogen uptakes in carbon nanostructures has been 
reported [93]. 

Fullerenes are synthesized carbon molecules that resemble a football in shape. 
Experiments showing adsorption and desorption of hydrogen (up to 6wt%) have been 
reported. The cyclic behavior of these materials is not impressive at this stage (as noted in 
[86]). 

Two types of carbon nanotubes have been explored: single-wall nanotubes (SWNT) 
and multi-wall nanotubes. Hydrogen can be stored in carbon nanotubes by chemisorption 
or physisorption mechanisms. In chemisorption the hydrogen is covalently bonded 
with host materials. The desorption temperatures are too high (> 500K). Physisorption 
occurs due to Van Der Waal’s attractive forces between hydrogen and the host material 
(carbon) that stabilizes the hydrogen on the surface of carbon nanotubes. Desorption 
temperatures of hydrogen for SWNTs, in the vicinity of room temperature, have been 
reported [94]. 

Activated carbons demonstrate high surface areas and are able to adsorb hydrogen in 
bulk material micro-pores [86]. Experimental data reported by numerous researchers have 
been tabulated for hydrogen storage (wt%) in carbon nanotubes, carbon nanofibers— 
graphitic planes that are intercalated at a fixed distance (0.355) from each other—and 
activated carbons [93]. 

The effective integration of carbon nanostructure hydrogen storage with RHGUS or 
automobiles will require the development of this material with a capability of operating 
near ambient conditions (25< T< 100°C and P> 1 atm). Effective modeling tools must be 
developed to predict accurately the mechanisms of hydrogen storage in carbon 
nanostructures. The cyclic performance of carbon nanostructures needs to be addressed 
along with the degradation of performance with time. 


2.8. Residential application 


Energy consumption of a house can be generally categorized in two forms, electrical 
and thermal (as a fuel). The energy demand of any particular house is a direct function of 
the type of materials of construction and their ability to transfer heat from/to the 
surrounding atmosphere. It is crucial to have an energy efficient design for future 
residential houses because it will not only reduce the overall energy demand of the house 
but also be helpful to design economically feasible RHGUS configurations for this 
application. 

The best possible scenario for designing RHGS/RHGUS for this application is to use 
actual energy demand data for the house on an annual basis. If this is not available, then 
the hourly electrical energy demand for the house can be obtained as given in what follows. 
Firstly, the cooling load and heating load demand can be estimated, as a function of design 
conditions and type of materials used, for the particular building [48]. The data for a 
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variety of structures and materials and their R-values (wall structure, roof, basement, 
doors, windows, lights, number of occupants, etc.) can be found using techniques 
summarized by ASHRAE [95]. Secondly, an overall heat transfer coefficient for the house, 
based on floor area, can be estimated for the design conditions of a particular location. 
This overall heat transfer coefficient can be further used to calculate the hourly electrical 
energy demand of the house as a function of ambient temperature. As an example, hourly 
energy demand for a standard house (1600 ft*) located in Las Vegas, NV, has been 
reported [48]. 

Once the hourly electrical energy demand is known, the next step is to obtain the 
thermal energy demand. The hydrogen-based car as an example can currently go 
approximately an average of 56.5miles/kg of hydrogen gas [96]. Assuming a certain 
percent of home refueling (0-100, where, 0 denotes the case with no hydrogen generated at 
home is used for fueling the car and 100 denotes all fuel needs of car are satisfied by 
hydrogen generated at home) can give the hydrogen required for driving the car. Besides 
this, the gas requirements for water heating, cooking, and drying can be estimated from 
data given by EERE [97]. Hourly residential gas requirements can be obtained by dividing 
the total amount of gas required by the number of hours in a year. 

Recently, the potential for research to design energy efficient buildings has been 
reported. Numerous aspects of energy efficient design of walls, windows, lighting, etc. have 
been discussed [98]. Another study reported the evaluation of evaporative-water cooled 
air-conditioning systems as compared with conventional air-cooled systems for a house in 
a desert climate [99]. 

Based on these electrical and gas requirements a RHGS/RHGUS system can be 
designed and its performance can be predicted for a particular house for any location. 
Economic issues for a house system designed according to the information given above are 
not addressed here. However, it is worth mentioning that if the house is already grid 
connected, then the grid can be easily used as short-term storage as compared with a 
battery. Certainly, this would not be a case for an off-grid application. 


3. Concluding assessment of state of system understanding 


From this review, it can be concluded that a large amount of effort has been taken in the 
recent past to model the various components of RHGS/RHGUS, either based on 
individual components or collectively for the overall system. Because of this vast literature 
on simulation and experimental efforts the question that might be asked is “Do we have a 
good understanding of the methodology to model the individual components as well as the 
overall system associated with a RHGS/RHGUS?” This article is an attempt to explore 
details of an answer to this question. Although simulation information presented in this 
paper is chosen to reflect the proper physics involved with the major components, further 
research is needed to reduce/eliminate generally the empirical nature of the models. 

The input energy resources needed for modeling purposes are well-quantified and over 30 
years’ of data are readily available for many locations in US. PV modeling discussed here is 
well developed and can be applied to the different types of cells currently being 
manufactured. Electrolyzer models are less well developed. For example, the ones 
applicable to PEM technology are currently in crude form and need attention to include 
more of the nature of the physics involved. Fuel cell modeling has received several years of 
careful attention, and this area is fairly mature. Recently published papers illustrate that 
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general performance aspects of these devices can be predicted quite well. The physics for 
the description of hydrogen storage of compressed gas are well understood. Only simple 
simulation tools have been put forth for describing metal hydride storage; more data are 
needed on large-scale applications to model accurately this method of containing the gas. 
Carbon-based storage technology is in very early development stages so its ability to store 
effectively on large scale is still in question. On the other hand, residential utilization of 
hydrogen can be readily represented in simulations. The energy requirements for a house 
can be modeled easily using methodology that is well developed. However, the current 
designs of buildings that, in general, require large amounts of energy, will have to have 
large RHGUS installations, and these will be very expensive. More energy efficient design 
of buildings will open this market more quickly. 


4. Recommendations for future work 


Based on the conclusions, some of the key aspects from a modeling standpoint for the 
RHGS/RHGUS that needs to be addressed by the scientific community are listed below: 


e Application of RHGS/RHGUS to homes is limited both by the relative immaturity of 
the technology, as well as the relatively large systems that would be needed. The 
development of standard designs of energy efficient residential houses is needed to 
reduce their energy demand. This will assist in decreasing the sizes of energy conversion 
equipment substantially. Appropriate methods to include the energy (electrical and 
thermal) efficient equipment and their effect on overall energy demands needs to be 
quantified. 

e Simulation approaches to represent the performance of photovoltaic cells are well 
developed. Quantifying the degradation of performance with time for these devices 
needs to be addressed. 

e Modeling insights for electrolyzers have to be improved significantly to show clearly the 
activation, ohmic and concentration polarization losses inside the cell. Since there is not 
enough data available for PEM electrolyzers, these aspects certainly should be 
addressed. Models need to be developed so that actual mechanisms of electrolysis can be 
reflected in it instead of theoretical estimates. 

e CFD modeling has been widely used to simulate the multi-physical phenomena 
associated with fuel cells. Not nearly as much attention has been given to the PEM 
electrolyzers, and the understanding of this technology could gain a great deal from 
these types of efforts. 

e Work related to the hydrogen storage technologies has been reported in [85]. The 
physics of some aspects of this important component of renewable systems are well 
understood, but there are several newer approaches that may hold great promise that 
need to be investigated further. 

è The degradation of performance validated with actual experimental data needs to be 
addressed for almost all the components of RHGUS. 

e Since the technology, manufacturing, and marketing of many of the RHGUS 
components are in early stages of development, costs of all the components of the 
RHGUS are very uncertain. This topic should be addressed now with prices estimated 
for the situation when this technology becomes more mature, and then revisited 
periodically as the market continues to grow. 
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